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1 Introduction 

The standard model (SM) of electroweak (EW) interactions HHSl posits the existence of the 
Higgs boson, a scalar particle associated with the field responsible for spontaneous EW sym¬ 
metry breaking I3H9I. The mass of the boson is not predicted by the theory In 2012 the ATLAS 
and CMS collaborations at the CERN LHC reported the observation of a new boson with a 
mass of about 125 GeV HUHTH. Throughout this paper, we denote the observed Higgs boson 
as h(125). Subsequent studies of the production and decay rates [[TsHSSI arid of the spin-parity 
quantum numbers of the new boson show that its properties are compatible with 

those expected for the SM Higgs boson. 

The observation of a Higgs boson with a mass of 125 GeV is also consistent with the unitarity 
constraints on diboson scattering at high energies ll28Vl37l . Nevertheless, there is a possibility 
that the newly discovered particle is part of a larger Higgs boson sector and thus only partially 
responsible for EW symmetry breaking. This can be realized in several scenarios, such as two- 
Higgs-doublet models IMl |3^ , or models in which the SM Higgs boson mixes with a heavy 
EW singlet Il40] - l63l . which predict the existence of additional resonances at high mass, with 
couplings similar to those of the SM Higgs boson. 

Previous searches at the LHG for heavy SM-like Higgs bosons have been reported by ATLAS 
and GMS. Based on a dataset of 4.7fb^^ at \/s = 7TeV and up to 5.9 fb^^ at ^/s = 8TeV, and 
combining all channels listed in Ref. |101, ATLAS excludes a SM-like heavy Higgs boson of 
131 < wjh < 559 GeV at 95% GL lUOl . The GMS collaboration reported a search in the WW and 
ZZ decay channels using an initial dataset of 5.1 fb~^ at \/s = 7 TeV and 5.3 fb^^ at \/s = 8 TeV, 
searching in the mass range 145 < nin < 1000 GeV, and excluding Higgs boson masses up 
to 710 GeV at 95% GL UMl . In this paper, we report on an extension of this search using the 
complete 7 and 8 TeV dataset. In addition, the search is interpreted in the context of the SM 
expanded by an additional EW singlet. Both the SM-like heavy Higgs boson as well as the EW 
singlet are denoted as H here. 

The analysis uses proton-proton collision data recorded with the GMS detector, corresponding 
to integrated luminosities of up to 5.1 fb^^ at \/s = 7TeV and up to 19.7fb^^ at \/s = 8TeV. 
The analysis is performed in a mass range 145 < nin < 1000 GeV exploiting both the H —)■ WW 
and H —> ZZ decay channels, with the lower boundary being chosen to reduce contamination 
from h(125). In the case of a Higgs boson decaying into a pair of W bosons, the fully leptonic 
(H —> WW —)• £v£v) and semileptonic (H —> WW —> £vqq) final states are considered in the 
analysis. Eor a Higgs boson decaying into two Z bosons, final states containing four charged 
leptons (H ^ ZZ ^ 2£2f), two charged leptons and two quarks (H ^ ZZ —> 2£2q), and two 
charged leptons and two neutrinos (H —> ZZ —> 2£2v) are considered, where £ = e or fi and 
£' = e, ji, or T. 

2 CMS detector 

The central feature of the GMS apparatus is a superconducting solenoid of 6 m internal di¬ 
ameter, providing a magnetic field of 3.8 T. Within the solenoid volume are a silicon pixel 
and strip tracker, a lead tungstate crystal electromagnetic calorimeter (EGAL), and a brass and 
scintillator hadron calorimeter, each composed of a barrel and two endcap sections. Muons 
are measured in gas-ionization detectors embedded in the steel flux-return yoke outside the 
solenoid. Extensive forward calorimetry complements the coverage provided by the barrel 
and endcap detectors. The first level of the GMS trigger system, composed of custom hardware 
processors, uses information from the calorimeters and the muon detectors to select the most 
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interesting events in a fixed time interval of less than 4 jis. The high level trigger processor farm 
further decreases the event rate from around 100 kHz to less than 1 kHz, before data storage. A 
more detailed description of the detector as well as the definition of the coordinate system and 
relevant kinematic variables can be found in Ref. 1651 . 


3 Signal model and simulations 

Several Monte Carlo event generators are used to simulate the signal and background event 
samples. The Higgs boson signal samples from gluon fusion (ggF, gg —> H), and vector bo¬ 
son fusion (VBF, qq ^ qqH), are generated with POWHEG 1.0 Il66ti70l at next-to-leading order 
(NLO) and a dedicated program HTTI used for angular correlations. Samples of WH, ZH, and 
tfH events are generated using the leading-order (LO) pythia 6.4 Il72l program. At the gener¬ 
ator level, events are weighted according to the total cross section P‘(pp —> H) 1731. which con¬ 
tains contributions from ggF computed to next-to-next-to-leading order (NNLO) and next-to- 
next-to-leading-log (NNLL), and from weak-boson fusion computed at NNLO. The WW(ZZ) 
invariant mass, mww i^zz), lineshape is affected by interference between signal and SM back¬ 
ground processes. The simulated OTh lineshape is therefore corrected to match theoretical pre¬ 
dictions Il74ti761 using the complex-pole scheme for the Higgs boson propagator. The procedure 
for including lineshape corrections and uncertainties from interference of the signal with back¬ 
ground processes for both ggF and VBF production are described below in the discussion of 
the lineshape corrections applied for the EW singlet interpretation. 

The background contribution from qq —)■ WW production is generated using MadGraph 
5.1 1771 , and the subdominant gg —> WW process is generated at LO with GG2ww 3.1 ITSlI . The 
qq —> ZZ production process is simulated at NLO with POWHEG, and the gg —> ZZ process is 
simulated at LO using GG2zz 3.1 173 . Other diboson processes (WZ, Z-|-jets, and 

W-|-jets are generated with PYTHIA and MadGraph. The tf and tW events are generated at 
NLO with POWHEG. For all samples, PYTHIA is used for parton showering, hadronization, and 
underlying event simulation. For LO generators, the default set of parton distribution func¬ 
tions (PDF) used to produce these samples is GTEQ6L |80l, while GTIO ISTl is used for NLO 
generators. The tau lepton decays are simulated with TAUOLA 1 82 1 . The detector response is 
simulated using a detailed description of the GMS detector, based on the Geant4 package 1831 . 
with event reconstruction performed identically to that of recorded data. The simulated sam¬ 
ples include the effect of multiple pp interactions per bunch crossing (pileup). The PYTHIA 
parameters for the underlying events and pileup interactions are set to the Z2 (Z2*) tune f84l 
for the 7 (8) TeV simulated data sample, with the pileup multiplicity distribution matching the 
one observed in data. 

The data are analysed to search for both a beyond the standard model (BSM) case in the form 
of an EW singlet scalar mixed with the recently discovered Higgs boson, h(125), as well as a 
heavy Higgs boson with SM-like couplings. The couplings of the two gauge eigenstates (h(125) 
and EW singlet) are phenomenologically constrained by unitarity and the coupling strength of 
the light Higgs boson is therefore reduced with respect to the SM case. The unitarity constraint 
is ensured by enforcing -|- C ^ = 1, where C and C' are defined as the scale factors of the 
couplings with respect to the SM of the low- and high-mass Higgs boson, respectively. The EW 
singlet production cross section is also modified by a factor /i' and the modified width is F'; 
they are defined as 


fi’ = C'2 (1 - Bne^), 


(1) 



where jBnew is the branching fraction of the EW singlet to non-SM decay modes. An upper limit 
at 95% CL can be set indirectly as C ^ < 0.28 using the signal strength fits to the h(125) boson 
as obtained in Ref. Il85l . 

This paper focuses on the case where C'^ < (1 — Bnew)- In this regime the new state is expected 
to have an equal or narrower width with respect to the SM case. Results are presented distin¬ 
guishing between the Bnew = 0 and Bnew > 0 cases. Under this hypothesis, signal samples 
with different Higgs boson widths are generated, scanning the C ^ and B^ew space. We follow 
the recommendations of the "LHC Higgs Cross Section Working Group" II73I described below. 

The SM signal mass lineshape generated with POWHEG is weighted in order to simulate the nar¬ 
row scalar singlet lineshape. Lor the ggL production mode, the weights are calculated using 
either the GGlZZ generator for the ZZ channel, or the POWHEG and MCFM 6.2 |86l generators 
for the interference calculation for the WW channel. Lor the VBL production mode, the in¬ 
terference weights are computed using the Phantom 1.2 USTlI generator, where the signal-only 
lineshape at LO is weighted based on results obtained with MadGraph generator predictions. 
The weights are defined as the ratio of the sum of a narrow resonance signal plus interference 
and the standard model signal lineshape as generated. The contribution from the interference 
term between the BSM Higgs boson and the background is furthermore assumed to scale ac¬ 
cording to the modified coupling of the Higgs boson as {}i + f)BSM = + kuC', where 

}i{I) is the signal strength (interference) in the BSM or SM cases. This assumption is based on 
the hypothesis that the couplings are similar to the SM case and simply rescale due to unitarity 
constraints. Systematic uncertainties considered for this procedure are detailed later. 

If the new resonance has a very small width, its production will tend to interfere less with 
the background continuum. Thus in the most interesting region of low-C'^, the effect of the 
interference and its exact modeling is of limited importance. Any possible interference between 
h(125) and its EW singlet partner |l88l|89l is assumed to be covered by a conservative systematic 
uncertainty. In addition to the EW singlet, the analysis searches for a heavy Higgs boson that 
gets produced and decays like the SM Higgs boson, but has a higher mass and interferes with 
h(125). 

4 Event reconstruction 

GMS uses a particle-flow (PL) reconstruction algorithm ll90ll9T] to provide an event description 
in the form of particle candidates, which are then used to build higher-level objects, such as 
jets and missing transverse energy, as well as lepton isolation quantities. Not all the channels 
considered here use the same selection criteria for their objects, but the common reconstruction 
methods are listed below. 

The high instantaneous luminosity delivered by the LHG provides an average of about 9 (21) 
pileup interactions per bunch crossing in 7 (8) TeV data, leading to events with several possi¬ 
ble primary vertices. The vertex with largest value of the sum of the square of the transverse 
momenta (px) for the associated tracks is chosen to be the reference vertex. According to sim¬ 
ulation, this requirement provides the correct primary vertex in more than 99% of both signal 
and background events. 

Muon candidates are reconstructed by using one of two algorithms: one in which tracks in the 
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silicon tracker are matched to hits in the muon detectors, and another in which a combined fit 
is performed to signals in both the silicon tracker and the muon system. Other identification 
criteria based on the number of measurements in the tracker and in the muon system, the 
fit quality of the muon track, and its consistency with its origin from the primary vertex are 
also imposed on the muon candidates to reduce the misidentification rate |92|. In some of 
the channels, vetos are imposed on additional low-px muons in the event, whose tracks in the 
silicon tracker fulfill more stringent requirements, but whose primary vertex association can be 
more relaxed. We call these muons soft muons. 

Electron candidates are reconstructed from superclusters, which are arrays of energy clusters 
along the cp direction in the ECAL, matched to tracks in the silicon tracker. A complemen¬ 
tary algorithm reconstructs electron candidates by extrapolating measurements in the inner¬ 
most tracker layers outward to the ECAL. Trajectories from both algorithms are reconstructed 
using the Gaussian sum filter algorithm 1931 , which accounts for the electron energy loss by 
bremsstrahlung. Additional requirements are applied to reject electrons originating from pho¬ 
ton conversions in the tracker material. Electron identification relies on a multivariate discrim¬ 
inant that combines observables sensitive to the bremsstrahlung along the electron trajectory, 
the geometrical and momentum-energy matching between the electron trajectory and the as¬ 
sociated supercluster, as well as ECAL shower shape observables. Electron candidates with 
a pseudorapidity tj of their ECAL supercluster in the transition region between ECAL barrel 
and endcap (1.4442 < \ri\ < 1.5660) are rejected, because the reconstruction of electrons in this 
region is compromised Il94l . 

The PE candidates are used to reconstruct hadronic tau candidates, Th, with the "hadron plus 
strip" algorithm Il95ll that is designed to optimize the performance of Th identification and re¬ 
construction by considering specific Th decay modes. The neutrinos produced in all t decays 
escape detection and are ignored in the Th reconstruction. The algorithm provides high Th iden¬ 
tification efficiency, approximately 50% for the range of Th energies relevant for this analysis, 
while keeping the misidentification rate for jets at the level of ~ 1%. 

Leptons produced in the decay of W and Z bosons are expected to be isolated from hadronic 
activity in the event. Isolation is defined as the scalar sum of the transverse momenta of the 
PE candidates found (excluding the selected leptons in the event themselves) in a cone of ra¬ 
dius R = 'J = 0.4, with (p being measured in radians, built around each lepton. 
We require that the isolation sum is smaller than 20% (15%) of the muon (electron) transverse 
momentum. To account for the contribution to the isolation sum from pileup interactions and 
the underlying event, a median energy density is determined on an event-by-event basis using 
the method described in Ref. Il96l . Eor electron candidates, an effective area that is proportional 
to the isolation cone, is derived to renormalize this density estimate to the number of pileup 
interactions, and is subtracted from the isolation sum. Eor muon and tau candidates, the cor¬ 
rection is done by subtracting half the sum of the transverse momenta of the charged particles 
not associated to the primary vertex in the cone of interest. Soft muons do not need to fulfill 
isolation requirements. 

The combined efficiency of lepton reconstruction, identification, and isolation is measured us¬ 
ing observed Z decays and ranges between 90% and 97% for muons, between 70% and 90% for 
electrons, and approximately 50% for hadronic taus depending on the pj and rj of the leptons. 

Jets are reconstructed from PE candidates 1971 using the anti-kj clustering algorithm |98| with 
a distance parameter of 0.5 (called AK5 jets), as implemented in the fastjet package 1991 , and 
with the Cambridge-Aachen algorithm 11001 with a distance parameter of 0.8 (called CAS jets); 
when not otherwise specified we use the term jets to mean AK5 jets throughout this paper. 
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Any reconstructed jet overlapping with isolated leptons within a distance of 0.5 (0.8) for AK5 
(CAS) jets is removed in order to avoid double counting the lepton as a jet. AK5 (CAS) jets are 
required to have | 7 /| < 4.7 (2.4). At hadron level, the jet momentum is defined as the vectorial 
sum of all particle momenta in the jet, and is found in the simulation to be within 5-10% of the 
true momentum over the whole pj range and detector acceptance. A correction is applied to 
the jet Pj to take into account the extra energy clustered in jets due to additional proton-proton 
interactions within the same bunch crossing ll^llOll . Other jet energy scale (JES) corrections 
applied are derived from the simulation, and are calibrated with in situ measurements of the 
energy balance of dijet and Z/ 7 +jets events. For some of the channels in this analysis, a combi¬ 
natorial background arises from low-px jets from pileup interactions, which get clustered into 
high-px jets. At v^s = 8 TeV the number of pileup events is larger than at = 7 TeV, and a mul¬ 
tivariate selection is applied to separate jets stemming from the primary interaction and those 
reconstructed due to energy deposits associated with pileup interactions II 102 I . The discrimi¬ 
nation is based on the differences in the jet shapes, on the relative multiplicity of charged and 
neutral components, and on the different fraction of transverse momentum, which is carried 
by the hardest components. Within the tracker acceptance, the tracks belonging to each jet are 
also required to be compatible with the primary vertex. 

At high Higgs boson masses, the pj of the Z or W boson is high enough that the two quarks 
from the vector boson decay are expected to be reconstructed as a single CAS jet, or merged 
jet. To improve background rejection and jet mass resolution, we apply a jet pruning algo¬ 
rithm I1103II . Additionally, the "N-subjettiness ratio" variable T 2 /Ti, a measure of the compat¬ 
ibility of a jet having N = 2 subjets, is used to reduce the backgrounds Ill03[ll04t . We require 
Tz/ti < 0.5 to avoid contamination from jets originating from the hadronization of gluons and 
single quarks. 

In some of the channels included in this analysis, the identification of jets originating from b 
quarks is important. These b jets are tagged with dedicated algorithms II105I , which are applied 
either directly to the AK5 jets, or to the subjets of the merged jets. The b-tagging algorithm 
used is either the Combined Secondary Vertex algorithm. Jet Probability algorithm, or Track 
Counting High Efficiency algorithm depending on the channel. Tagged b-jet candidates are 
required to have pj > 30GeV and to be within the tracker acceptance (|?/| < 2.4). For jets in 
this kinematic range with a b-tagging efficiency of 70%, the misidentification probability from 
light quark and gluon jets is approximately 1 %. 

The missing transverse momentum vector, is defined as the projection on the plane per¬ 

pendicular to the beams of the negative vector sum of the momenta of all reconstructed parti¬ 
cles in an event 1191111061 . Its magnitude is referred to as 


5 Data analysis 

The results presented in this paper are obtained by combining searches exploiting different 
Higgs boson production and decay modes as detailed in Table All final states are exclusive, 
without overlap between channels. For the WW —> £v£v and ZZ —> 2£2£' channels a mass 
range starting from 110 GeV has been analyzed in other contexts (e.g. HU), but in this paper 
both analyses are restricted to searches above 145 GeV. In the rest of this section, the individ¬ 
ual analysis strategies and details are defined including a discussion of the leading systematic 
uncertainties. In the next section, statistical interpretations of the individual and combined 
searches are given for the Higgs boson with SM-like couplings hypothesis as well as the EW 
singlet extension of the SM. 
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Table 1: A summary of the analyses included in this paper. The column "H production" indi¬ 
cates the production mechanism targeted by an analysis; it does not imply 100% purity in the 
selected sample. The main contribution in the untagged categories is always ggP. The (jj)vBF 
refers to a dijet pair consistent with the VBF topology, (jj)w(z) ^nd (J)w(z) refer to a dijet pair 
and single merged jet from a Lorentz-boosted W (Z) with an invariant mass consistent with a 
W (Z) dijet decay, respectively. The superscript "0,1,2 b tags" refers to the three possible cat¬ 
egories of b tag multiplicities. Exclusive final states are selected according to the lepton and 
reconstructed jet content of the event. The mass range under investigation and the mass reso¬ 
lution are also listed. Mass ranges differ with the sensitivities of each channel. 


H 

H 

Exclusive 

No. of 

mu range 

mu 

decay mode 

production 

final states 

channels 

[GeV] 

resolution 

WW ^ iv(.v 

untagged 

((ee, pp), ep) + (0 or 1 jets) 

4 

145-1000|f| 

20% 


VBF tag 

((ee, pp), ep) + (jj)vBF 

2 

145-1000^ 

20% 

WW ^ ^vqq 

untagged 

(ev, pi/) + (jj)w 

2 

180-600 

5-15% 


untagged 

(ev, pv) + (J)w + (0+1-jets) 

2 

600-1000 

5-15% 


VBF tag 

(ev, pv) + (J)w + (ij)vBF 

1 

600-1000 

5-15% 

zz ^ im' 

untagged 

4e, 4p, 2e2p 

3 

145-1000 

1-2% 


VBF tag 

(4e, 4p, 2e2p) + (jj)vBF 

3 

145-1000 

1-2% 


untagged 

(ee, pp) + (Th+h, TeTh, T^^Th, T^Ty) 

8 

200-1000 

10-15% 

ZZ ^ U2v 

untagged 

(ee, pp) + (0 or > 1 jets) 

4 

200-1000 

7% 


VBF tag 

(ee, pp) + (jj)vBF 

2 

200-1000 

7% 

zz ^ 2£2q 

untagged 

(ee,pp) + (jj)"T 2 btags 

6 

230-1000 1 

3% 


untagged 

(ee, pp) + (J)”’'/"-*- 

6 

230-1000" 

3% 


VBF tag 

(ee, pp)+ (jj)z^'^'’*'‘®'+ (jj)vBF 

6 

230-1000" 

3% 


VBF tag 

(ee, pp) + (J)7;^'^^‘^s% (jj)vBF 

6 

230-1000" 

3% 


“EW singlet model interpretation starts at 200 GeV to avoid contamination from h(125). 
*'600-1000 GeV for = 8 TeV only. 

"For -^s = 8 TeV only. 


5.1 H ^ WW ^ iviv 

In the H —>• WW —)■ iv(,v channel the Higgs boson decays to two W bosons, both of which decay 
leptonically, resulting in a signature with two isolated, oppositely charged, high-px leptons 
(muons or electrons) and large due to the undetected neutrinos. A complete description 
of the analysis strategy is given in Ref. IITSl . For this analysis, we require triggers with either 
one or two high-px muons or electrons. The single muon or electron triggers are based on 
relatively tight lepton identification with px thresholds from 17 to 25 GeV (17 to 27 GeV) in the 
muon (electron) channel. The higher thresholds are used for periods of higher instantaneous 
luminosity. The dilepton trigger px thresholds for the leading and trailing leptons were 17 and 
8 GeV, respectively. 

Gandidate events must contain two reconstructed leptons with opposite charge, px > 20 GeV 
for the leading lepton, and px > 10 GeV for the subleading one. Only muons (electrons) with 
Ip I < 2.4 (2.5) are considered in this channel. The analysis is very similar to that in the Higgs 
boson discovery lTlJ[l2l, but additionally uses an improved Higgs boson lineshape model. 

Events are classified into three mutually exclusive categories, according to the number of re¬ 
constructed jets with Px > 30 GeV. The categories are characterized by different signal yields 
and signal-to-background ratios. In the following, these are referred to as 0-jet, 1-jet, and 2-jet 
multiplicity categories. Events with more than two jets are considered only if they are con¬ 
sistent with the VBF signature: a dijet mass of the two highest px jets greater than 500 GeV, 
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I > 3.5, and no additional jets with pj > 30 GeV in the p region between these two leading 
jets. Signal candidates are further divided into same-flavor (SF) dilepton or e^e^) and 

different-flavor (DF) dilepton (p^e^) categories. The bulk of the signal arises through direct W 
boson decays to muons or electrons, with the small contribution from W —> rv —^ X decays 
implicitly included. 

In addition to high-px isolated leptons and minimal jet activity, significant is expected 
to be present in signal events, as opposed to none to moderate in background. For this 
channel, an variable is employed, defined as (i) the magnitude of the component 

transverse to the closest lepton, if A(p(i, £™'®®) < 7 r/ 2 , or (ii) the magnitude of £^^®® otherwise. 
This observable more efficiently distinguishes Z/ j* —)■ t+t“ background events in which the 
F^i®® is preferentially aligned with the leptons, and Z/ 7 * —)■ events with mismeasured 

£^^®®. Since the £^p® resolution is degraded as pileup increases, the minimum of two different 
observables is used: the first includes all particle candidates in the event, while the second 
uses only the charged particle candidates associated with the primary vertex. Events with 
£^p® > 20 GeV are selected for this analysis. 

The backgrounds are suppressed using techniques described in Refs. 1TT1|121. Top quark back¬ 
ground is controlled with a selection based on the presence of a soft muon and b-jet tag¬ 
ging II105I . Rejection of events with a third lepton passing the same requirements as the two 
selected leptons reduces both WZ and W 7 * backgrounds. 

The Drell-Yan (DY) process produces SF lepton pairs (p^p^ and e+e^) and therefore addi¬ 
tional requirements are applied for the SF final state. Firstly, the resonant component of the 
DY background is rejected by requiring a dilepton mass outside a 30 GeV window centered 
on the Z boson mass. The remaining off-peak contribution is further suppressed by requir¬ 
ing significant missing transverse energy in the event. For the = 7TeV data, we require 
£^p® > (37 -|- Nvtx/2) GeV directly, while for = 8 TeV data, a boosted decision tree (BDT) 
multivariate discriminant including £^p® is used. For events with two jets, the dominant source 
of misreconstructed £“'®® is the mismeasurement of the hadronic recoil, and optimal perfor¬ 
mance is obtained by requiring £^^®® > 45 GeV. Finally, the momenta of the dilepton system 
and of the most energetic jet must not be back-to-back in the transverse plane. For the \/s = 
7 TeV data, we require the angle to be less than 165 degrees , while for ^ys = 8 TeV data, this in¬ 
formation is included in the BDT. These selections reduce the DY background by three orders 
of magnitude, while rejecting less than 50% of the signal. 

The final analysis in the DF and in the SF 0-jet or 1-jet categories is performed using a two- 
dimensional fit in two observables, the dilepton invariant mass m££ and the transverse mass 
nij determined between the transverse dilepton system, pj, and the £^^®®: 

mj = y2p"£“i®®(l - cos AO(pf, £^'®®)), (3) 

where £^'®®) is the azimuthal angle between the dilepton transverse momentum and 

the £^®®. 

Figure shows the distributions of these two observables for the 0-jet DF category. The fit 
ranges in ni££ and nij are dependent on the nin hypothesis. For niu > 200 GeV, additional 
selections requiring pj > 50 GeV for the leading lepton and rtij > 80 GeV are imposed to 
suppress the h(125) contribution. A cross-check counting analysis is performed by applying 
selection criteria to several kinematic observables including ni££ and nij. 


For the VBF production mode 11107111 lOl . the cross section is roughly ten times smaller than for 
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Figure 1: Distributions of mu (left) and mj (right) for the 0-jet DF category of the H —)■ WW —> 
£viv search. The uncertainty in the background histograms includes the systematic uncertain¬ 
ties on all background estimates and is centered on the sum of all backgrounds, including the 
h(125) in red. The W -|-jets distributions include the contributions from QCD multijet processes 
as well. The red open histogram shows five times the expectation for a = 400 GeV SM-like 
Higgs boson. The selection has been optimized to suppress the h(125) contribution as explained 
in the text. 

ggF at lower hypotheses and is roughly three times smaller at the highest wjh hypothesis. 
We optimize the selection in the 2-jet category to tag these VBF-type events by requiring the 
mass of the dijet system to fulfill wjjj > 500 GeV, and the angular separation of the two jets to 
pass |A//jj| > 3.5. Given the small event yield in this category, the signal extraction in the DF 
category is performed using a one-dimensional fit mmu where an dependent requirement 
on the transverse mass is imposed. A counting analysis is performed in the SF category and is 
used as a cross-check in the DF category 

The normalization of the background contributions relies on observed events rather than simu¬ 
lation whenever possible, and exploits a combination of techniques IfTTlIT^ . The tt background 
is estimated by extrapolation from the observed number of events with the b-tagging require¬ 
ment inverted. The DY background measurement is based on extrapolation from the observed 
number of and e+e^ events with the Z boson veto requirement inverted. The background 

of W -|- jets and Quantum Ghromodynamics (QGD) multijet events is estimated by measuring 
the number of events with one loosely isolated lepton. The probability for such loosely isolated 
nongenuine leptons to pass the tight isolation criteria is measured in observed data using mul¬ 
tijet events. In the 0-jet and 1-jet bins, the nonresonant WW contribution is estimated from a fit 
to the data while in the 2-jet bin it is estimated from simulation. Other backgrounds, such as 
V+Z/ 7 * and triple boson production (WV) are estimated from simulation and are small. 

Experimental effects, theoretical predictions, and the choice of event generators (powheg, 
GG2ww, MadGraph, Phantom) are considered as sources of systematic uncertainty, and 
their impact on the signal efficiency is assessed. The overall signal normalization uncertainty 
is estimated to be about 20%, and is dominated by the theoretical uncertainty associated with 
missing higher-order QGD corrections and PDF uncertainties, estimated following the PDF4LHG 
recommendations IIHTl IlllblMl . The total uncertainty in the background estimation in the 
H —> WW signal region is about 15% and is dominated by the statistical uncertainty in the 
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observed number of events in the background control regions. 

5.2 H ^ WW ^ ivqq 
5.2.1 Unmerged-jet category 

In the H —> WW —)■ £vqq channel we search for a Higgs boson decaying to WW, where one 
W decays leptonically, thus providing a trigger handle for the event, while the other decays 
hadronically This channel has a larger branching fraction than the two-lepton final state and 
allows one to reconstruct the Higgs boson candidate invariant mass II115I . The final state con¬ 
sists of an isolated electron or muon, and two separated jets. The main experimental 

challenge is to control the large W+jets background. 

We use data collected with a suite of single-lepton triggers mostly using pj thresholds of 24 
(27) GeV for muons (electrons). Basic kinematic selections are applied to the final-state objects 
to reduce the background contribution. The muon (electron) is required to be isolated and have 
Pt > 25 (35) GeV. A veto is imposed on additional muons and electrons in the event to reduce 
backgrounds from DY events. The events are required to have > 25 (30) GeV in case 

of muons (electrons). The transverse mass nij = V2pjE^^^^{l — cos AO), with AO being the 
azimuthal angle between the lepton pj and needs to be greater than 30 GeV. The two 

leading-px jets in the event each must have px > 30 GeV and must together have an invariant 
mass, mjj, between 66 and 98 GeV. An unbirmed maximum likelihood fit is performed on mjj to 
determine the background normalizations in the signal region for each Higgs mass hypothesis 
independently Events that contain a b-tagged jet are vetoed in order to reduce the background 
from top quark decays. 

To further exploit the differences in kinematics between the signal and background, a likeli¬ 
hood discriminant is constructed using angles between the Higgs boson decay products that 
fully describe the Higgs boson production kinematics IIT9l ItTII and the lepton charge. The lep¬ 
ton charge provides discrimination power because of the charge asymmetry in the background, 
which is not present in the signal. This approach improves the expected sensitivity to a Higgs 
boson across the entire mass range. 

The background normalization in the signal region is extracted for each Higgs boson mass hy¬ 
pothesis independently with an unbinned maximum likelihood fit to the dijet invariant mass 
distribution, mjj, of the two leading jets. The signal region corresponding to the W boson mass 
window, 66 < wjjj < 98 GeV, is excluded from the fit. The background is overwhelmingly due 
to events from W+jets production. The normalization of the W+jets component is a free param¬ 
eter and is determined in the fit. The electroweak diboson, tt, and single top shapes are based 
on simulation and their normalizations are constrained to the theoretical predictions with as¬ 
sociated uncertainties. In the case of diboson backgrounds, both WW and WZ normalizations 
come from NLO predictions. The uncertainties in the normalization of the W+jets component 
obtained from the fit, as well as those from all other backgrounds, are included in the limit 
calculation as systematic uncertainties. 

For the signal search, we use the binned distribution of /«ww/ which is computed using a neu¬ 
trino longitudinal momentum, pz, that is determined from the constraint that the leptonically 
decaying W boson is produced on-shell. The ambiguity of two possible solutions is resolved 
by taking the solution that yields the smaller \pz\ value for the neutrino. According to sim¬ 
ulation over 85% of signal events are assigned the correct pz value, thus improving the mass 
resolution, especially at low Higgs mass. The contribution of each background in the signal 
region (66 < my < 98 GeV) is described by an expected shape and normalization. The distri- 
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bution of the W+jets contribution is parameterized with a polynomial functional form deter¬ 
mined through simulation. The parameters of the function are determined from the observed 
data spectrum in the binned likelihood fit, which is also used to determine the exclusion lim¬ 
its. The other background shape distributions are parameterized using simulation, and their 
normalizations and uncertainties are included as nuisance parameters. Figure depicts the 
njww distribution for the muon category for two different Higgs boson mass hypotheses. The 
observed data are compared to background and signal expectations. 
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Figure 2: The WW invariant mass distribution with the fit projections in the signal region 
66 < mjj < 98GeV, for the muon channel of the unmerged-jet category. The V+jets back¬ 
ground includes both, the large contribution from W+jets production, and the much smaller 
component of Z+jets. The blue curve on the left (right) shows 50 (5) times the expectation for a 
mil = 200 (500) GeV SM-like Higgs boson. 


Experimental effects, theoretical predictions, and uncertainties due to the choice of fit functions 
are considered as sources of systematic uncertainty. Because of the large background, the dom¬ 
inant source of systematic uncertainty is the shape uncertainty of the W+jets mww distribution, 
followed by the normalization uncertainty that is extracted from the my fit. The main uncer¬ 
tainty in the signal normalization stems from the uncertainty in the efficiency of the likelihood 
discriminant selection. This effect occurs because of mismodeling of the likelihood discrimi¬ 
nant and is studied in a signal-depleted region of the analysis. 


5.2.2 Merged-jet category 

In the highest mass range of this search, from 600 to 1000 GeV, the pj of the decaying W bosons 
is typically greater than 200 GeV. At this pi, the daughter quarks of the hadronically decay¬ 
ing W are often merged into a single jet to the point where traditional dijet searches cannot 
be performed. For a signal mass of 600 GeV (1 TeV), and signal events falling in the detector 
acceptance, approximately 65% (82%) of the hadronic W decay products are contained in a 
cone of AR < 0.8; alternatively, approximately 10% (42%) of the hadronic W decay products 
are separated by a distance of AR < 0.5 and would not be reconstructed by the standard GMS 
AK5 jet finding algorithm. The larger cone GA8 jet affords more signal acceptance in the single 
jet signature while not losing events when the decay products are separated by AR < 0.5. In 
this case, we use jet substructure techniques for identifying jets that have originated from a 
hadronically decaying W boson with high Lorentz boost. 
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As in the unmerged case, we use data collected with the single electron or muon trigger. The 
dominant background is W+jets with a smaller background contribution from tf. Remaining 
backgrounds arise from WW, WZ, ZZ, and single top quark production. 

The hadronic W boson candidate is reconstructed with CAS jets to increase acceptance. To 
reduce the contribution of quark- and gluon-initiated jets from QCD processes, a selection is 
made on the pruned jet mass of the CAS jet of 65 < wjj < 105 GeV and the N-subjettiness ratio 
Tz/ti < 0.5. 

The kinematic selections of muons and electrons are slightly more stringent than those for the 
unmerged case, with the pj thresholds of isolated muons and electrons being 30 and 35 GeV, 
respectively. The requirement is increased to be greater than 50 (70) GeV in the muon 
(electron) channel. The pj of both, leptonically and hadronically decaying W boson candidates, 
is required to be greater than 200 GeV in order to select events with a large Lorentz boost. 
Additional selections are made to ensure that the W boson candidates are sufficiently separated 
in a back-to-back topology As in the unmerged category, there is an additional requirement to 
veto events with a b tagged jet in order to reduce the background from top quark decays. 

To increase the sensitivity of the analysis, the analysis is split into exclusive jet multiplicity 
categories: 0+1-jet and 2-jets. These additional AK5 jets must pass a pj threshold of 30 GeV and 
be separated from the hadronic W candidate in the event by a distance AR > 0.8. Additionally, 
the 2-jet category has further requirements to better identify a topology consistent with the VBF 
production mode. The two highest pj AK5 jets in the event are required to satisfy | A/^jj | > 2.5 
and > 250 GeV. There is also a requirement on the invariant mass of the W boson candidates 
and the nearest AK5 jet to be greater than 200 GeV in order to reject other top quark background. 
The final discriminating distribution is the wJww shape reconstructed from the lepton + 

+ GAS jet system. Because the number of events in the 2-jet category is limited, the muon and 
electron datasets are merged. 

The background is estimated using observed data rather than simulation wherever possible. 
The dominant W+jets background normalization is determined from a fit to the pruned jet 
mass sideband where other backgrounds contribute less significantly The W+jets shape is 
determined from extrapolating the pruned jet mass sideband region into the signal region. The 
tf background is estimated by inverting the veto on AK5 b-tagged jets, which yields a high- 
purity tf control sample. The normalization of tf in the signal region is then corrected by a 
data-to-MG scale factor determined from the tf control sample. Finally the tf control sample 
is also a good source of W bosons with high Lorentz boost, which are used to calibrate the jet 
substructure selection efficiency 

Figure displays the final ?«vvw distributions after all selections and background estimations 
for the 0+1-jet category for the muon channel only on the left and for the 2-jet category on the 
right. In the 2-jet category, there is an excess observed around 750 GeV with local significances 
of 2.6 a at 700 GeV and 2.1 a at 800 GeV. No such excess is observed in the 0+1-jet multiplicity 
bin. 

Sources of systematic uncertainty are similar to the unmerged-jet category The dominant ex¬ 
perimental uncertainties are the background normalization and shape uncertainty, particularly 
in the 2-jet category, where the number of events is small, as well as the uncertainty in the signal 
cross section due to the hadronic W candidate selection. The dominant theoretical uncertainty 
comes from the uncertainty in the cross section when dividing the analysis into exclusive jet 
multiplicity bins. 
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^WW [G®V] ^WW [G®V] 

Figure 3: The final WW invariant mass distribution is shown for the 0+1-jet bin category for the 
muon channel only (left) and for the 2-jet bin category (right). Points represent the observed 
data, shaded graphs represent the background and dashed graphs represent five times the 
expectation for a wjh = 800 GeV SM-like Higgs boson from ggF and VBF production, separately 

5.3 H ^ ZZ ^ lilt 

This analysis seeks to identify Higgs boson decays to a pair of Z bosons, with one Z decaying 
to a pair of muons or electrons (Z —> 21, with £ = or e), and the second decaying to electrons, 
muons or taus (Z —)■ It, with t = e,}i or t) IfTTlIT^ITT^ . This channel has extremely low back¬ 
ground, and the presence of four leptons in the final state allows reconstruction and isolation 
requirements to be very loose, leading to a high selection efficiency This channel is one of the 
most sensitive channels across the entire mass range. 

For this analysis, we require triggers with two high-px muons or electrons. The dilepton trig¬ 
ger px thresholds for the leading and trailing leptons were 17 and 8 GeV, respectively. Events 
included in the analysis contain Z boson candidates formed from a pair of leptons of the same 
flavor and opposite charge. Decay muons or electrons are required to be isolated, and to orig¬ 
inate from the primary vertex. Muons (electrons) are required to have px > 5(7) GeV and 
Ip I < 2.1 (2.5), while taus are required to have a visible transverse momentum px > 20 GeV 
and IpI < 2.3. We reconstruct the Z —> tt in the following decay modes: Z —> ThTh, Z ^ TeTh, 
Z ^ T^T^, and Z —> TeT^,. Overlap with the channel is avoided by excluding events with 
both taus decaying to electrons or muons. 

For the 2121 final state, the lepton pair with invariant mass closest to the nominal Z boson 
mass, denoted Zi, is identified and retained if 40 < niz-^ < 120 GeV. The second Z boson, 
denoted Z2, is then constructed from the remaining leptons in the event, and is required to 
satisfy 12 < mz 2 < 120 GeV. If more than one Z2 candidate remains, the ambiguity is resolved 
by choosing the leptons of highest px. Amongst the four candidate decay leptons, at least one 
should have px > 20 GeV, and another should have px > 10 GeV. This requirement ensures 
that selected events correspond to the high-efficiency plateau of the trigger. 

For the 2£2t final state, events are required to have one Zi —)■ 21 candidate, with one lepton 
having px > 20 GeV and the other px > 10 GeV. The leptons from leptonic decays of the tau 
are required to have px > 10 GeV. The invariant mass of the reconstructed Zj candidate is 
required to satisfy 60 < < 120 GeV. The Z2 candidate mass reconstructed from the visible 
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tau decay products (visible mass), m-c^, must satisfy 30 < < 90GeV for events with at 

least one hadronically decaying tau, and 20 < nirr < 90 GeV for events with two leptonically 
decaying taus. Events with both taus decaying to muons or electrons are excluded in order to 
avoid overlap with the 2i2£ channel. The selections on the reconstructed Zi mass for the 2£2t 
final state are tighter than for the 2£2£ charmel, because the search range starts at ?«h = 200 GeV 
rather than = 145 GeV. 

In order to further separate signal from background and to distinguish different signal produc¬ 
tion mechanisms, we use a matrix element likelihood approach, which relies on probabilities 
for an event to come either from signal or background using a set of observables, such as an¬ 
gular and mass information, which fully characterize the event topology in its center-of-mass 
frame |711 . 

The events are categorized according to their jet multiplicity, counting jets with pj > 30 GeV. 
In the 0- or 1-jet category the pj spectrum of the four-lepton system is exploited to distinguish 
between ggF and vector boson induced production modes, such as VBF and associated produc¬ 
tion with a vector boson (VH). The events selected with two or more jets potentially come from 
several sources: background (predominantly ZZ+2 jets), VBF signal, VH signal (with V—>-2 
jets), and ggF signal gg —> H-|-2 jets. In order to isolate the production mechanism, we use 
a matrix-element-based probability fl/Tl 11171 for the kinematics of the two jets and the Higgs 
boson to come either from the VBF process or from the gg —> H-|-2 jets process. This discrimi¬ 
nant is equally efficient to separate VBF from either H-|-2 jets signal or ZZ+2 jets background, 
because both processes show very similar jet kinematics. 

To allow estimation of the tf, Z+jets, and WZ+jets backgrounds, a Zi+X control region is de¬ 
fined, well separated from the expected signal region. In addition, a sample of Z-i + £ events, 
with at least one reconstructed lepton in addition to a Z, is defined in order to estimate the 
lepton misidentification probability, which is the probability for non-prompt leptons and other 
particles, which are not leptons, to be reconstructed as leptons and to pass the isolation and 
identification requirements. The contamination of the sample from WZ events containing a 
genuine additional lepton is suppressed by requiring the energy imbalance measured in the 
transverse plane to be below 25(20) GeV for the 2£2£ (2£2r) channel. The event rates mea¬ 
sured in the background control region are extrapolated to the signal region. The systematic 
uncertainties associated with the reducible background estimate vary from 30% to 100% and 
are combined in quadrature with the statistical uncertainties. 

The cross section for ZZ production at NFO is calculated with mcfm. The theoretical uncer¬ 
tainty in the cross section is evaluated as a function of varying both the renormalization 

and factorization scales and the PDF set, following the PDF4FHG recommendations IIHTl IlllT 

[m. 

The reconstructed invariant mass distribution for 2£2£, combining the 4/^, 2}i2e, and 4/^ chan¬ 
nels, is shown in Fig. (left), compared with the expectation from SM background processes. 
Figure]^ (right) plots the reconstructed visible mass distribution for the 2£2t selection, combin¬ 
ing all 2£2t final states. 

The background shapes are taken from simulation, with rates normalized to the observed data. 
The measured Z boson and ZZ masses are underestimated due to the undetected neutrinos 
in tau decays. To compensate for this effect we correct the momenta of the visible tau decay 
products to constrain the Z boson mass to 91.2 GeV. The correction is applied on the selected 
events and only affects the final mass shape of the ZZ system. The observed mass distributions 
are well-matched to the SM background expectation. 
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m^z [GeV] [GeV] 

Figure 4: Distribution of the four-lepton reconstructed mass for the sum of the 4fi, 2}i2e, and 
4e channels (left), and for the sum over all 2i2T channels (right). Points represent the observed 
data, shaded histograms represent the background, and the red open histogram shows the 
expectation for a wjh = 350 GeV SM-like Higgs boson. 

In the 2£2l channel limits on the production of heavy Higgs bosons are extracted using the 
unbrnned four lepton mass distribution and the correlation of the kinematic discriminant with 
the Higgs boson mass. In the case of the tagged 2-jet category a third dimension is introduced 
using the correlation of the aforementioned VBF discriminant with the four-lepton mass. In the 
untagged 0/1-jet category the transverse momentum of the four lepton system is used in place 
of the VBF discriminant. For the 212t final state, limits are set using the mun distribution. 

5.4 H ^ ZZ ^ 2ilv 

This analysis seeks to identify Higgs boson decays to a pair of Z bosons, with one Z boson 
decaying to neutrinos and the other decaying to leptons. The analysis strategy is based on 
selection requirements in the nij) phase space, with selections adjusted for different mu 

hypotheses I1118L Here, the transverse mass is determined between the transverse dilepton 
system, pj, and the 

As in the previous Section, we use data collected with the trigger requiring two high-px elec¬ 
trons or muons. Events are required to have a pair of well-identified, isolated leptons of same 
flavor (p^p^ or e^e~), each lepton with px > 20 GeV, with an invariant mass within a 30 GeV 
window centered on the Z boson mass. The px of the dilepton system is required to be greater 
than 55 GeV. The presence of large (70 GeV or more, depending on wjh) in the event is 
also required. 

To suppress the Z+jets background, events are rejected if the angle in the transverse plane 
between the and the closest jet with px > 30 GeV is smaller than 0.5 radians. Events where 

the lepton is mismeasured are rejected if > 60 GeV and A(p{£, £™^®®) < 0.2. The top quark 
background is suppressed by applying a veto on events having a b tagged jet with px > 30 GeV 
and IpI < 2.5. To further suppress this background, a veto is applied on events containing an 
additional, soft muon with px > 3 GeV, typically produced in the leptonic decay of a b quark. 
To reduce the WZ background, in which both bosons decay leptonically, any event with a third 
lepton (p or e) with px > 10 GeV, and passing the identification and isolation requirements, is 
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rejected. 

The search is carried out in two mutually exclusive categories. The VBF category contains 
events with two or more jets in the forward region, with a > 4 requirement between 

the two leading jets, and with the invariant mass of those two jets greater than 500 GeV. In 
addition, the two leptons forming the Z candidate are required to lie between these two jets in 
rj, while no other selected jets with px > 30 GeV are allowed in this region. The ggF category 
includes all events failing the VBF selection, and is subdivided into subsamples according to 
the presence or absence of reconstructed jets with pj > 30 GeV. The event categories are chosen 
in order to maximize the expected cross section limit. In the case of the VBF category, a constant 
£miss ^ 70 QgY requirement are used, as no gain in sensitivity is obtained with a 

selection dependent on the Higgs boson mass hypothesis. In the case of ggF, we apply an 
?«H-dependent lower limit on ranging from 80 to 100 GeV over the search range. 

The background composition is expected to vary with the ?«h hypothesis. At Iow-wjh, Z+jets 
and tf processes are the largest contributions, while at wjh > 400 GeV the irreducible ZZ and 
WZ backgrounds dominate. The ZZ and WZ backgrounds are estimated from simulation and 
are normalized to their respective NLO cross sections. The Z+jets background is modeled from 
a control sample of events with a single photon produced in association with jets. This proce¬ 
dure yields a good model of the distribution in Z+jets events. 

The uncertainty associated with the Z+jets background is affected by residual contamination 
of the 7 +jets control sample from processes involving a photon and genuine £^^®®. We do not 
explicitly subtract this contamination, but include a shape uncertainty in the final fit, which is 
allowed to vary this small residual background between 0 and 100 % of this estimate. 

Background processes that do not involve a Z boson (nonresonant background), include tf, 
single top quark production, W+jets and WW, and are estimated with a control sample of 
DF dilepton events (p^e^) that pass the full event selection. This method cannot distinguish 
between the nonresonant background and a possible contribution from H —> WW —> 2£2v 
events, which are treated as part of the nonresonant background estimate. The nonresonant 
background in the or e+e^ final states is estimated by applying a scale factor to the 

selected f<^e^ events, obtained from the events in the sidebands of the Z boson peak (40 < 
mi£ < 70 GeV and 110 < m££ < 200 GeV). The uncertainty associated with the estimate of this 
background is determined to be 25%. No significant excess of events is observed over the SM 
background expectation. 

The mj and £^^®® distributions used to extract the results are shown in Fig.|^for all three event 
categories with the muon and electron channels combined. 

5.5 H ^ ZZ ^ 2£2q 

The H —)■ ZZ —)■ 2^2q channel has the largest branching fraction of all H —)■ ZZ charmels under 
consideration, but also a large background contribution from Z+jets production. Furthermore, 
the reconstruction of the hadronically decaying Z boson is more difficult than the fully leptonic 
final states. The dominant background is due to the DY process and can be reduced using 
b-tagging requirements on the selected jets, as will be described later in this section. 

As in the other channels with a Z boson in the final state, we use data collected with the trigger 
requiring two high-px muons or electrons. Reconstructed muons and electrons are required to 
have px > 40(20) GeV for the leading (second) lepton. Muons are required to have |p| < 2.4, 
and electrons |p| < 2.5. Jets are required to have px > 30GeV and |p| < 2.4. Each pair of 
oppositely charged leptons of the same flavor and each pair of jets are considered as Z boson 
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Figure 5: The final transverse mass (left, center) and missing transverse energy (right) distri¬ 
butions are shown for all three event categories of the H —> ZZ —> 2£2v channel: (left) events 
with zero jets, (center) events with at least one jet, but not passing the VBF selection, (right) 
VBF events. The expected distributions from the different background processes are stacked 
on top of each other. The red open histograms show the expectation for a ma = 400 GeV SM- 
like Higgs boson separating the ggF and VBF contributions. In the case of the 0-jet category, 
the VBF contribution is multiplied by a factor of 10 to increase visibility. 


candidates. 


To increase the sensitivity to a possible signal, the main analysis is complemented with dedi¬ 
cated selections for the VBF signature. VBF events are identified requiring two additional jets 
with px > 30 GeV, mjj > 500 GeV, and |Apjj| > 3.5. 


Analogously to the H —> WW —> £i/qq channel, the analysis is modified for very high Higgs 
boson masses (wjh > 600 GeV) to account for the fact that the two jets originate from a Lorentz- 
boosted Z boson and therefore they may be reconstructed as a single, merged jet (Sections 
and 5.2.2| . Information about the internal structure of this kind of jet II103II is used in or¬ 
der to gain some insight about the origin of the jet, distinguishing the DY background from 
jets produced from boosted Z bosons. To reduce the contamination from nonboosted Z -|- X 
backgrounds, the selection for the merged-jet topology requires the hadronically decaying Z 
boson to have pj > 100GeV and \rj\ < 2.4, and the leptonically decaying Z boson to have 
Pt > 200 GeV. 


In order to exploit the different jet composition of signal and background, events are classified 
into three mutually exclusive categories according to the number of selected b tagged jets: 0 
b tag, 1 b tag, and 2 b tag. This distinction is not done in the VBF oriented selection where 
the main discrimination is given by a multivariate discriminator for VBF topology based on 
angular and energy information of the two VBF tag jets. 

Background contributions are reduced by requiring 71 < Wjj < 111 GeV and 76 < < 

106 GeV in the selected events. The presence of Z bosons decaying to leptons and dijets makes 
this selection very efficient for signal, whereas the continuous background gets largely reduced. 
In case of the merged-jet analysis, the dijet requirement is applied on the merged-jet mass, after 
applying the pruning procedure. 

An angular likelihood discriminant is used to separate signal-like from background-like events 
in each category IItTI 11191 . In order to suppress the substantial expected tf background in the 
2 b tag category, a discriminant is used, defined as the logarithm of the likelihood ratio of the 
hypothesis that the is equal to the value measured by the PF algorithm and the null hy¬ 
pothesis = 0 lEl- This discriminant provides a measure of whether the event contains 
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genuine missing transverse energy. When an event contains multiple Z boson candidates pass¬ 
ing the selection requirements, those with jets in the highest b tag category are retained for 
analysis. If multiple candidates are still present, the ones with mjj and niu values closest to 
the Z boson mass are retained. In the case of the merged-jet category all the requirements on 
jets with respect to b tagging and the likelihood discriminant are applied to the two subjets 
reconstructed inside the merged jet. In case of events with a merged-jet Z boson candidate in 
addition to a dijet one, the candidate from the merged jet gets selected. 

The dominant Z+jets background is estimated using simulated events properly corrected to 
reproduce the yield of the observed data in the control regions, defined as 60 < ?«jj < 71 GeV 
and 111 < njjj < 130 GeV. Simulated events are weighted to reproduce the pj spectrum of 
the 11]] system in these control regions. In case of the merged-jet category both the pj of the 
and systems have weights applied. The normalization of this distribution is taken from 
observed data and used as an additional constraint for this background. 

Other backgrounds are estimated using simulated events. These include diboson and top quark 
events. In the VBF selection, the background from top quark events is small. In the ggF analysis 
it is estimated from observed data using a control sample of ep events, invoking lepton flavor 
symmetry that is satisfied by several background contributions but not in the signal events. 
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Figure 6: The nizz invariant mass distribution after final selection in three categories of the 
H ^ ZZ ^ 2£2q dijet channel: 0 b tag (left), 1 b tag (center), and 2 b tag (right). Points 
with error bars show distributions of observed data. Solid histograms are depicting the back¬ 
ground expectation from simulated events for the different components. The red open his¬ 
togram shows the expectation for a jmh — 400 GeV SM-like Higgs boson. 

The distributions of nizz in the signal region are shown in Fig. |^for the three b tag multiplicities 
of the dijet category comparing observed data with background expectations. Good agreement 
is observed within the uncertainties. For the merged and VBF categories good agreement is 
also observed. The dominant systematic uncertainties are due to the b-tagging performance 
and the JES Il97l . Further systematic effects are due to the uncertainty in the predicted signal 
and background shapes used in the analysis. The distributions of mzz in the signal region are 
used to extract the final signal and background yields. 

5.6 Systematic uncertainties 

Systematic uncertainties for the various final states come from common treatment of the signal 
model assumptions, reconstructed objects used in the analysis and a few common experimental 
effects. 

Uncertainties on the cross section for the production of heavy Higgs bosons arise from un¬ 
certainties in the combined choice of the PDFs and ocs, as well as in the renormalization and 
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5 Data analysis 


factorization scales Il73l , which are typically 6-7% and 7-12%, respectively, for the ggF produc¬ 
tion mechanism, and 1-2% and 2-5%, respectively, for production via VBF. Additionally, we 
add an uncertainty in the background coming from off-shell h(125) production, which we esti¬ 
mate with GG2zz (Phantom) for the ggF (VBF) case. We find that at the largest ntn values, the 
size of the effect is approximately 3% of the total background. Uncertainties on the signal line- 
shape reweightmg with interference varies for the ggF and VBF modes. For ggF, we follow the 
prescription in Ref. 1731 . which considers the NNLO contribution to the signal interfering with 
the gg —> ZZ background process. For VBF, without a full prescription, we assign systematic 
uncertainties coming from renormalization and factorization scale variations in the PHANTOM 
generator. 

Other common systematic effects come from the luminosity uncertainty, which is 2.2% (2.6%) 
for the 7 (8) TeV data. Uncertainties on the muon and electron reconstruction efficiencies, and 
JES and jet energy resolution (JER) are correlated among the various final states, where all these 
effects are subdominant. The lepton fake rate is largest for the H —> ZZ —> 2£2f channel, in 
which we consider fake leptons at relatively lower px than in the other charmels. A summary 
of the systematic uncertainties per channel is given in Table 

Table 2: Sources of systematic uncertainties considered in each of the channels included in this 
analysis. Uncertainties are given in percent. Most uncertainties are affecting the normalisation 
of the observed data or simulated events, but some are uncertainties on the shape of kinematic 
distributions. Wherever ranges of uncertainties are given, they are either ranges in mn, jet 
multiplicity categories, or dependent on the production mode. 


Source of uncertainty 

WW 

— >■ £v£v 

WW 

^^vjj 

WW 

-^£v] 

H ^ ZZ 
^ 2£2£' 

H ^ ZZ 
^ 2£2v 

ZZ 

^2f2q 

Experimental sources 


Luminosity, 7 (8) TeV 

2.2 (2.6) 

2.2 (2.6) 

2.2 (2.6) 

2.2 (2.6) 

2.2 (2.6) 

2.2 (2.6) 

i trigger, reco, id, iso 

1-4 

1-2 

1-2 

0.5-7 

2-3 

1.8-2 

i mom./energy scale 

2-4 



0.5-30 

1-2 

0.1-0.4 

£ misid. rate 




30 



JES, JER, 

2-35 

<1 

2 

5-30 

1 

1-13 

Pileup 


<1 



1-3 

1 

b-tag/ mistag 



2.5 


1-3 

1-6 

W-tag/Z-tag 



7.5 



0-9.3 

Signal selection eff. 


10 

2 




Monte Carlo statistics 

1-20 




1-2 

0-6 

Background estimates 


tf,tw 

20 

7 

6-30 


25 
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Z+jets 
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100 

16 

ZZ 
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12 


W+jets 

40 
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8 


25 


WW 

8-30 

10 

30 


25 


WZ, Wy* 

3-50 


30 
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Theoretical sources 


^(gg ^ H) 

10-13 
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11-13 

10-13 

10-13 

10-13 

u(qq ^ H) 

2.6-5.8 

2.6-3.6 

3.6-5.8 

2.6-5.8 

2.6-5.8 

2.6-5.8 

H lineshape 




5 

2-8 

0-7 

H-WW (ZZ) interference 

1-27 
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Jet binning 

7-35 


7-35 
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6 Statistical interpretation 

The combination of the measurements in the different channels presented in this paper requires 
the simultaneous analysis of the data selected by all individual analyses, accounting for all sta¬ 
tistical and systematic uncertainties, as well as their correlations. The statistical methodology 
used in this combination was developed by the ATLAS and CMS Collaborations in the context 
of the LHC Higgs Combination Group lH^ 11201 Il21l . Upper limits on the model parameters 
are set for different wjh hypotheses using a modified frequentist method referred to as the CLg 
method II1221I1231 , where a likelihood ratio test statistic is used in which the nuisance param¬ 
eters are profiled. In the likelihood ratio, the total number of observed events is compared to 
the signal and background predictions by means of a product of Poisson probabilities. The 
predictions are subject to the multiple uncertainties described in the previous section, and are 
handled by introducing nuisance parameters with probability density functions. The nuisance 
parameters modify parametrically the expectations for both signal and background processes. 
Furthermore, a signal strength modifier (ji) is used to scale the Higgs boson cross sections of all 
production mechanisms by the same factor with respect to their SM predictions while keeping 
the decay branching fractions unchanged. 

6.1 SM-like Higgs boson search 

The combined results obtained for a heavy Higgs boson with SM-like couplings for all the 
different contributing final states are displayed in Fig. On the left, the observed 95% CL limit 
is shown for each final state. The expected combined 95% CL limit of the six channels is plotted 
as a dashed black line, while the yellow shaded region is the ±2u uncertainty in the expected 
limit. On the right, the expected and observed limits are displayed for each of the individual 
channels as well as the combined result. The top right panel shows the WW final states, while 
the ZZ final states are displayed in the bottom right panel. In the lower mass region of the 
search range, the most sensitive channels are H —> ZZ —)■ AE and H —)■ WW —> 2E2v. At the 
highest masses, the H —> ZZ —> 2£2v channel has the best sensitivity, while H —)■ ZZ —)■ AE, 
H —)• WW —)• 2E2v, H —> WW —> £i/qq, and H —)■ ZZ —>• 2£2q contribute significantly. 

Features in the combined observed limit can be traced to corresponding features in the limits of 
the individual channels. At lower masses below 400 GeV, there are oscillations in the observed 
limit due to the high resolution channel, H —)■ ZZ —> AE, and the narrow width of the heavy 
Higgs boson in this mass range. An excess in the combined limit at around 280 GeV is related 
to a small excess in the channels H —> ZZ —>• AE and H —> WW —> 2E2v. The small excess 
of observed events seen around 700 GeV in the H —)■ WW EvJ merged-jet category is not 
supported by the other channels and is reduced to less than 0.5 cr in the combination. The com¬ 
bined upper limits at the 95% GL on the product of the cross section and branching fractions 
exclude a Higgs boson with SM-like couplings in the full search range 145 < wjh < 1000 GeV. 

6.2 EW singlet Higgs boson search 

We interpret the search results in terms of a heavy Higgs boson in the EW singlet extension 
of the SM. The parameters of the model are C'^, the heavy Higgs boson contribution to elec- 
troweak symmetry breaking, and Bnew, the contribution to the Higgs boson width of non-SM 
decays, and are defined in Section]^ Figure|^shows the expected and observed upper limits at 
95% GL on the singlet scalar cross section with respect to its expected cross section as a function 
of its mass. 

The region above the curves shows the parameter space that is expected to be excluded at 
95% GL. We show the exclusion region for various values of Bnew We find a large region 
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Figure 7: Upper limits at the 95% CL for each of the contributing final states and their combi¬ 
nation. The theoretical cross section, crgM, is computed in Ref. Il73l . The observed and expected 
limits of the six individual channels are compared with each other and with the combined re¬ 
sults (right), for H —)■ WW channels (top right panel) and H —> ZZ channels (bottom right 
panel) separately. 

of versus mass parameters to be excluded for various values of Bnew We also plot the 
Fh(i 25 ) = 1 ± 0.14 118511 indirect constraint C'^ < 0.28 at 95% CL for 13new=0. The upper dash- 
dotted line shows the cutoff of the allowed region for Snew = 0.5 where the width of the heavy 
Higgs boson becomes larger than the SM width at that mass hypothesis. 

In order to understand the constraints of these results in a model-independent approach, we 
further subdivide the results into categories. In Fig. we show the limits in various configu¬ 
rations. At the top of Fig. are the limits we obtain when we combine the ZZ (top left) and 
WW (top right) charmels separately. Since the ZZ channels are more sensitive in the search for 
a Higgs boson with SM-like couplings, they better constrain the BSM case as well. The bottom 
of Fig. 1^ shows the combined 95% CL for all final states but only the ggF or VBF production 
mechanism for the heavy Higgs boson. In the heavy Higgs boson with SM-like couplings sce¬ 
nario, we assume the ratio of the cross sections for various production mechanisms to be the 
same as in the SM case. 











6.2 EW singlet Higgs boson search 
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Figure 8: Upper limits at the 95% CL on the EW singlet extension. Upper limits are displayed 
as a function of the heavy Higgs boson mass and the model parameter C'^ for different values 
of i5new- The upper dash-dotted line indicates where, for Snew = 0.5, the variable width of the 
heavy Higgs boson reaches the width of a SM-like Higgs boson. The lower dash-dotted line 
displays the indirect limit at 95% CL on C'^ from the measurement of h(125). 
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CMS up to 5.1 fb'’ {7 TeV) + up to 19.7 fb'’ (8 TeV) 






Figure 9: Upper limits at the 95% CL on the EW singlet extension. Upper limits are displayed 
as a function of the heavy Higgs boson mass and the model parameter C'^ for different values 
of Bnew All considered ZZ decay charmels combined (top left). All considered WW decay 
channels combined (top right). Limits for the ggF production mode only (bottom left). Limits 
for VBF production only (bottom right). The dash-dotted line in the two left plots indicates 
where, for Snew = 0.5, the variable width of the heavy Higgs boson reaches the width of a 
SM-like Higgs boson. 






























23 


7 Summary 

Combined results are presented from searches for a heavy Higgs boson in H —> WW and H —> 
ZZ decay channels, for Higgs boson mass hypotheses in the range 145 < nin < 1000 GeV. In the 
case of a Higgs boson decaying into a pair of W bosons, the fully leptonic (H —> WW —)■ Ivlv) 
and semileptonic (H —> WW —> £vqq) final states are considered in the analysis. For a Higgs 
boson decaying into two Z bosons, final states containing four charged leptons (H —> ZZ —> 
2i2f), two charged leptons and two quarks (H —> ZZ —> 2£2q), and two charged leptons and 
two neutrinos (H —> ZZ —> 2£2v) are considered, where £ = eor and £' = e, }i, or t. 

The observed data are interpreted both in the context of a heavy Higgs boson with SM-like cou¬ 
plings and decays, as well as a search for a heavy, narrow resonance as an EW singlet partner 
of the SM Higgs boson at 125 GeV. No significant excess over the expected SM background has 
been observed and exclusion limits have been set. In the case of the search for a heavy Higgs 
boson with SM-like couplings and decays, we exclude the existence of such a heavy Higgs bo¬ 
son over the entire search range of 145 < OTh < 1000 GeV. For the EW singlet partner of the SM 
Higgs, the parameters of the model are C'^, the heavy Higgs boson contribution to EW symme¬ 
try breaking, and the contribution to the Higgs boson width of non-SM decays. We find 
that a large part of the C'^ versus mass parameter space is excluded for various values of B^ew 
Additionally, we present limits for the EW singlet model for different production mechanisms, 
ggE and VBE, and WW and ZZ decay modes separately. 
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